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ABSTRACT 


The  data  of  dye  patch  experiments  made  off  New  Jersey  and  on  the  Bahama 
Banks  are  analyzed.  The  second  moment  of  the  concentration  distribution  in¬ 
creases  with  square  of  time  in  the  initial  stage  and  with  cube  of  time  later, 
in  agreement  with  a  theory  of  a  two-particle  problem.  In  order  to  explain 
the  curved  trailing  part  of  the  dye  patches,  refraction  of  a  submerged  object 
is  treated  theoretically.  It  is  found  that  the  effect  of  refraction  may 
cause  a  curved  image,  but  the  curvature  due  to  refraction  is  much  smaller 
than  actually  observed.  The  vertically  differential  advection  as  the  main 
cause  of  the  curved  pattern  is  treated  analytically,  and  the  result  is 
applied  to  computing  vertical  eddy  viscosity  from  the  observed  curved 
trailing  portion  of  dye  patches.  The  mechanism  of  striated  patterns  of  dye 
patches  is  discussed  from  the  point  of  view  of  hydrodynamic  instability  and 
of  rolling  up  of  streak  lines. 


I.  INTRODUCTION 


As  emphasized  by  Frenkiel  (1958),  turbulence  is  a  normal  condition  of 
fluid  motion  encountered  in  nature;  and  laminar  motion,  to  which  more  than 
of  Sir  Horace  Lamb’s  exhaustive  Hydrodynamics  was  devoted,  is  rather 
exceptional.  Almost  any  water  movement  in  the  ocean,  from  the  general  circu¬ 
lation  to  sound  waves,  is  more  or  less  influenced  by  turbulence.  Therefore, 
it  is  rather  surprising  that  veiy  few  studies  have  been  done  on  oceanic  tur¬ 
bulence,  although  there  has  been  substantial  progress  in  recent  years,  both 
in  theoretical  and  experimental  researches  on  turbulence  in  general  and  to 
some  degree  atmospheric  turbulence  in  particular. 

Among  others,  there  are  two  main  reasons  for  the  lack  of  progress  in 
studying  oceanic  turbulence.  One  is  that  measurement  of  quantities  in  the 
ocean  pertinent  to  turbulence  is  complicated,  owing  to  technical  difficulties 
of  instrumentation  and  the  rigors  of  oceanographic  field  work.  The  other 
is  a  theoretical  difficulty.  A  theory  of  isotropic  homogeneous  turbulence 
was  well  developed  (Batchelor,  1953),  and  it  was  found  that  most  of  the 
results  are  valid  in  the  atmospheric  turbulence.  However,  in  the  ocean  such 
validity  is  dubious,  partly  due  to  the  horizontal  boundaries  and  free  surface 
both  of  which  have  no  counterpart  in  the  atmosphere,  and  partly  due  to  the 
condition  that  the  energy  of  turbulence  is  mostly  supplied  through  the  free 
surface  in  the  form  of  wind  stresses. 

However,  from  another  point  of  view,  these  liabilities  for  studying 
oceanic  turbulence  become  assets.  In  the  experimental  phase,  measurement  of 
turbulence  characteristics  by  Eulerian  methods  is  far  more  difficult  in  the 
ocean  than  in  the  atmosphere,  since  current  meters  capable  of  recording  small 
fluctuations  of  currents  in  the  sea  are  yet  to  be  made.  However,  Lagrangian 
methods  can  be  more  easily  utilized  in  the  ocean  than  in  the  atmosphere.  For 
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instance,  tracing  of  current  drogues  can  be  done  relatively  easily  for  a 
long  period  of  time  and  will  give  fairly  accurate  pictures  of  flow  patterns 
and  their  fluctuations,  while  tracking  of  pilot  balloons  in  the  atmosphere  is 
feasible  only  for  a  short  time  and  is  not  accurate  enough  to  determine  fluctu¬ 
ations  of  the  wind.  Diffusion  of  dye  patches  has  also  become  a  powerful  tool 
in  studying  horizontal  turbulence  in  the  ocean,  owing  to  recently  developed 
techniques  for  measuring  concentration  and  for  navigation.  Dye  patches 
can  be  followed  for  a  fairly  long  time  and  movements  of  patches  also  occur, 
in  most  cases,  almost  horizontally  with  a  small  spreading  along  the  vertical 
direction*  These  are  remarkable  advantages,  considering  the  case  of  diffusion 
of  smoke,  since  one  cannot  make  experiments  of  large  scale  diffusion  by  using 
smoke  more  than  a  few  hours  without  facing  a  charge  of  willful  contamination 
of  the  air* 

As  for  complications  arising  from  the  existence  of  horizontal  boundaries 
and  free  surface,  we  may  consider  that  the  ocean  provides  an  ideal  natural 
laboratoiy  for  turbulence  with  different  scales  from  several  meters  to  thou¬ 
sands  of  kilometers  (Frenkiel,  1962).  Deviation  from  the  homogeneous,  iso¬ 
tropic  turbulence  in  the  ocean  poses  challenging  problans  to  theoreticians. 

The  experiments  described  in  this  report  were  made  on  the  shelf  off .New 
Jersey  and  on  the  Bahama  Banks.  In  these  experiments  twenty  to  fifty  gallons 
of  solution  of  Rhodamine  B  dye  was  introduced.  The  concentration  of  dye  in 
the  patches  was  measured  continuously  with  a  standard  Turner  Model  III 
fluorometer,  while  a  research  vessel  made  several  crossings  of  the  patch 
within  a  period  of  a  few  hours.  Aerial  photographs  of  dye  patches  were 
taken  from  a  light  plane  at  ten  minute  intervals.  Technical  details  of 
dye  sampling  methods,  navigation  methods,  and  aerial  photography  have  already 
been  explained  in  a  previous  report  (Costin,  Davis,  Gerard  and  Katz,  1963), 
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which  is  referred  to  hereafter  as  ENY. 

This  paper  presents  analysis  and  theoretical  considerations  on  the 
experiments  in  two  areas,  on  the  Bahama  Banks  and  off  New  Jersey,  These  two 
areas  are  representative  of  oceanic  regions  with  different  turbulence  charac¬ 
teristics,  The  Bahama  Banks  are  extremely  shallow  with  uniform  depth,  and 
the  level  of  background  turbulence  (horizontal  components)  is  correspondingly 
low.  The  sea  off  New  Jersey  is  typical  of  the  Continental: Shelf,  and  the  back 
ground  turbulence  energy  is  higher  than  on  the  Bahama  Bank. 

II.  TWO  PARTICLE  MODEL  OF  DYE  PATCH  DIFFUSION 

The  results  of  diffusion  experiments  off  New  Jersey  (ENY)  and  on  the 
Bahama  Bank  (Costin,  1963)  indicated  that  the  patterns  of  dye  patches  were 
complicated  and  dissimilar  to  those  predicted  by  existing  diffusion  theories. 
Particularly,  striated  patterns  and  elongation  of  patches  are  notable  features 
which  are  not  compatible  with  a  model  of  diffusion  by  isotropic  turbulence. 
However,  when  the  logarithm  of  concentration  is  plotted  against  the  area 
enclosed  by  the  contour  of  that  concentration,  there  is  a  linear  relationship 
in  most  cases  of  these  experiments  (ENY).  Therefore,  the  procedure  of  using 
areas  enclosed  by  contours,  instead  of  linear  distances,  may  provide  some 
kind  of  averaging  process,  on  which  the  conventional  theories  of  turbulent 
diffusion  are  based.  When  equivalent  radii  r0  is  defined  as  equal  to  /a77l  > 
where  A  is  the  area  enclosed  by  a  contour,  the  relationship  between  concen¬ 
tration  C  and  re  is  expressed  by; 


CD 


p 

where  Cm  is  the  maximum  concentration  and  is  the  second  moment  of  the 


concentration  distribution. 
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Diffusion  from  an  instantaneous  source  which  is  moving  with  a  general 
current  is  defined  as  relative  diffusion  and  is  discussed  approximately  as 
two  particle  problem  (Batchelor,  19^0 ;  Gifford,  193>7)*  However,  it  is  neces¬ 
sary  to  discuss  this  type  of  diffusion  as  a  many  particle  problem.  The  diffi¬ 
culty  in  the  many  particle  problem  lies  in  the  fact  that  the  Lagrangian  corre¬ 
lation  functions  for  turbulent  velocities  are  not  yet  obtained.  A  simple 
reason  for  this  difficulty  can  be  explained  as  follows: 

Consider  that  N  particles  are  released  from  a  point  initially. 

It  is  assumed  that  there  is  no  mean  current.  The  x-coordinate  of  the 
i-th  particle  at  time  t  is  expressed  with  x^(t).  The  average  with  N 
particles  is  represented  by  the  bar.  Then  at  time  t. 


or: 


where  u-p  is  the  x-component  of  the  velocity  of  the  i-th  particle.  If 
all  the  particles  are  constrained  on  the  x-axis,  as  in  the  case  of  the 
ordinary  one  particle  problem,  (A  (  tY  l{  ( ~t  j  is  the  correlation  function 
which  depends  (t  -  )  only  in  homogeneous  turbulent  flow.  However,  the 

velocities  u^(t)  and  )  are  not  only  dependent  on  the  time,  but 

also  on  the  positions  of  the  i-th  particle  x^(t)  and  x^(^).  Therefore 
we  have  to  write: 
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Equation  (3)  indicates  that  the  correlation  function  u(t)u(/^')  is, 
indeed,  dependent  on  x^(t)  and  x^(^). 

Batchelor  (195>0)  treated  the  two  particle  problem  from  a  basis  of  dimen¬ 
sional  analysis  instead  of  going  into  details  of  the  Lagrangian  correlation 
(3).  According  to  his  theory,  the  relative  motion  of  two  particles  which 
are  released  at  the  distance  £*  Q  is  governed  by  turbulence  of  the  inertial 
subrange  if  the  Reynolds1  number  of  the  turbulence  is  sufficiently  large,  as 
in  the  ocean.  In  the  inertial  subrange,  energy  of  turbulence  is  transferred 
from  larger  eddies  to  smaller  ones  due  to  inertial  forces,  and  energy  of  the 
smaller  eddies  is  finally  dissipated  by  viscosity  at  the  rate  0  •  Therefore, 

physical  quantities  pertinent  to  the  process  in  this  range  are  Q  and  <00 , 

1  c 

since  viscosity  affects  only  the  smallest  eddies.  Then  dimensional  reasoning 


leads  to  the  relation: 


at 


=  Eta  H  ('  t  i) 
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»  1 

for  the  rate  of  change  of  mean  square  separation  £  (t)  ,  where  F  is  a  uni- 
versal  function*  When  t  is  sufficiently  small*  the  relative  velocity  of  the 
two  particles  remains  constant,  and  thus  the  left-hand  side  of  (It)  is  pro¬ 


portional  to  t  (Batchelor,  195*6),  so  that: 
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On  the  other  hand,  if  t  is  so  large  that  the  relative  motion  of  the  two 

particles  no  longer  depends  on  the  initial  separation  C  ,  equation  (It) 
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becomes: 
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Since  the  constant ^  ^  eqUation  (l)  is  proportional  to  p  ^(t) 
(Gifford,  195>7)*  validity  of  relations  (6)  and  (8)  can  be  tested  by  plotting 
determined  from  the  curves  of  concentration  versus  area,  against  time. 
The  result  is  shown  in  Figure  1,  in  which  the  values  of  ()F  determined  from 
the  experiments  off  New  Jersey  and  on  the  Bahama  Bank  are  plotted  in  logarith¬ 
mic  scale.  Equation  (6)  can  be  written  as: 

=  Co*  st.  (  t/'-t-  t  o  ) 


&2 


(t) 


(6') 


where  tQ  is  the  constant  which  is  proportional  to  0  3  ^  .  Since 

there  is  no  deductive  method  to  determine  constants  in  equations  (6)  and  (8), 
we  have  to  be  satisfied  with  comparing  the  trend  of  theoretical  curves  with 
the  values  determined  from  the  data.  Figure  1  indicates  that  the  rate  of 
increase  of ^  wj_th  time  agrees  with  the  theoretical  formula  (8)  for  a 
large  time  interval.  In  order  to  compare  equation  (61)  with  the  data,  curves 
are  plotted  for  t0  =  0,  200  min.,  300  min,  and  !?00  min.,  respectively. 

Since  t0  depends  on  the  initial  condition  of  dye  patches,  different  sets  of 
experimental  data  correspond  to  different  values  of  t0  in  the  first  stage, 
which  is  represented  by  equation  (6)  or  (61)*  It  can  be  seen  clearly  that 
there  is  a  second  stage  corresponding  to  (8),  except  in  the  data  of  the 
South  Bahama  Bank.  This  figure  also  indicates  that  the  turbulent  energy 
dissipation  Q  is  larger  in  the  area  off  New  Jersey  than  on  the  Bahama  Bank, 
since  the  values  of 0(  ^  are  larger  in  the  former  area  for  the  same  range  of 
t.  This  is  plausible,  because  the  shallow  water  of  the  Bahama  Bank  inhibits 
full  development  of  turbulence. 
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Since  the  constant ^  ^  in  equation  (l)  is  proportional  top  ^(t) 
(Gifford,  195>7),  validity  of  relations  (6)  and  (8)  can  be  tested  by  plotting 
determined  from  the  curves  of  concentration  versus  area,  against  time. 
The  result  is  shorn  in  Figure  1,  in  which  the  values  of  determined  from 
the  experiments  off  New  Jersey  and  on  the  Bahama  Bank  are  plotted  in  logarith¬ 
mic  scale*  Equation  (6)  can  be  written  as: 

fFTt)  =  co-v^t.  (e+ t*2-) 

where  tn  is  the  constant  which  is  proportional  to  ?  3  F  ^  .  Since 

there  is  no  deductive  method  to  determine  constants  in  equations  (6)  and  (8), 
we  have  to  be  satisfied  with  comparing  the  trend  of  theoretical  curves  with 
the  values  determined  from  the  data*  Figure  1  indicates  that  the  rate  of 
increase  of  q/  ^  with  time  agrees  with  the  theoretical  formula  (8)  for  a 
large  time  interval.  In  order  to  compare  equation  (6*)  with  the  data,  curves 
are  plotted  for  tQ  =  0,  200  min.,  300  min.  and  500  min.,  respectively. 

Since  t0  depends  on  the  initial  condition  of  dye  patches,  different  sets  of 
experimental  data  correspond  to  different  values  of  t0  in  the  first  stage, 
which  is  represented  by  equation  (6)  or  (6*)«  It  can  be  seen  clearly  that 
there  is  a  second  stage  corresponding  to  (8),  except  in  the  data  of  the 
South  Bahama  Bank.  This  figure  also  indicates  that  the  turbulent  energy 
dissipation  Q  is  larger  in  the  area  off  New  Jersey  than  on  the  Bahama  Bank, 
since  the  values  of 0(  ^  are  larger  in  the  former  area  for  the  same  range  of 
t.  This  is  plausible,  because  the  shallow  water  of  the  Bahama  Bank  inhibits 


full  development  of  turbulence. 
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III.  ELONGATION  AND  CURVATURE  OF  DIE  PATCHES 


In  many  cases  during  the  dye  experiments,  dye  patches  showed  a  tendency 
to  elongate  in  a  direction  of  surface  currents,  as  described  in  ENY.  Such 
elongation  may  be  explained  partly  by  anisotropy  of  turbulence,  and  partly 
by  vertically  differential  advection  of  the  mean  flow*  The  first  effect  was 
verified  experimentally  in  a  shear  flow  of  a  flume,  in  which  the  large  eddies 
are  elongated  in  the  direction  of  the  mean  flow  (Townsend,  19^6).  Since  the 
surface  current,  in  most  cases  of  the  experiments,  were  generated  by  surface 
winds,  there  was  a  vertical  shear  near  the  surface.  Thus,  turbulent  diffu¬ 
sion  was  usually  stronger  in  the  direction  of  the  surface  flow  than  in  the 
lateral  direction,  although  an  accurate  measurement  of  fluctuations  of  cur¬ 
rents  in  the  sea  is  yet  to  be  done  to  verify  this  effect. 

The  second  effect  was  usually  more  obvious  in  many  experiments.  Some 
aerial  photographs  clearly  indicate  that  the  leading  portion  of  a  patch  with 
the  highest  concentration  moves  faster  near  the  surface  while  the  trailing 
portion  lies  below  the  surface,  as  often  seen  from  a  band  of  dye  stirred  up 
by  crossing  of  a  research  vessel  (Example,  Figure  9  in  ENY).  Analytically, 
this  type  of  dye  distribution  is  expressed  by  a  solution  of  the  following 
diffusion  equation: 


where  S  is  the  concentration,  U(z)  andV(z)  are  horizontal  components  of  the 
mean  flow,  kx  and  Ay  are  horizontal  eddy  diffusivity,  and  the  vertical  eddy 
diffusion  is  neglected.  The  solution  of  equation  (l),  corresponding  to  the 
initial  source  concentration  F(x,y, z),  is  given  by: 
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For  the  line  source  at  the  origin  (  x  a  y  »  0)*  the  solution  is  given 
by  the  integrand  of  (2)  with  F  =  1  and  (\>»  L  =  0.  Further,  when  V(z)  a  0, 


o 


u. 


the  concentration  in  aerial  photographs  is  represented  by: 

h  (i,  %  z,t  )  e  d  z 

J0  2  »-j 

=(>Trt)''(A>A3)'2exp/'a  x< 


in  which  k  is  a  vertical  extenction  coefficient  of  dye  in  the  water,  and  H  is 
the  depth*  The  concentration  Sa  is  computed  from  equation  (3)  for  the  case 
of  constant  current  (U  a  U0)  and  for  the  case  of  vertically  variable  current 


t  X.  j  j  f  ■"  y\  x-y 

U0  t-o ->  by  taking  k  =  0.  The  result  is  schematically  shown 


in  Figure  2,  in  which  the  concentration  is  plotted  with  arbitrary  units.  The 
figure  indicates  that  the  dye  patch  becomes  elongated  owing  to  a  long  trailing 
part  after  a  certain  period. 


IF.  CURVATURE  CF  THE  PATCHES 

It  was  discovered  in  a  number  of  experiments  that  the  trailing  portion 
of  the  dye  patches  showed  a  pronounced  clockwise  curvature  from  the  source 
near  the  surface  to  the  tail  in  increasing  depth  when  the  wind  speed  was 
fifteen  knots  or  more  (Gerard  and  Katz,  1963 ,  hereafter  cited  as  GK). 

Gerard  and  Katz  considered  that  the  curved  tail  of  dye  patches  is  caused  by 
the  effect  of  an  Ekman  spiral  (actually,  effect  of  differential  advection  due 
to  the  Ekman  flow  in  the  upper  layer),  owing  to  the  fact  that  the  trail  is 
curved  clockwise  in  most  cases. 

Another  possible  explanation  for  this  phenomenon  is  the  effect  of  re¬ 
fraction  at  the  sea  surface  when  an  underwater  object  is  observed  from  above. 


FIGURE  2, 


Concentration  diagrams  for  an  instantaneous  point  source  either  in 
a  shear  flow  or  in  a  uniform  flow. 

Non-dimensional  parameters:  £  =  X  Vo  A_/ 

,  'f'  =  i  ijc  2 

The  groups  A,  B3  C,  D  and  E  represent  the  concentration  at  h*9,  9.7 

38*9,  5>8*3,  and  87«5>>  respectively. 


The  subscripts  1  and  2  of  each  group  indicate  the  concentration  for  the 
shear  flow  TJ  =  ~Uo  ^ 


rr  ^“7„  *  respectively. 


and  the  uniform  flow 
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in  order  to  determine  this  effect  on  curvature  of  the  image  at  the  sea 
surface,  the  geometrical  arrangement  of  the  object,  image  and  camera  (or  eye) 
is  shown  in  Figure  3*  In  this  figure  P(x,y,z)3  P*  (X,Y,0),  and  Q(x0,y0,z0)  are 
the  underwater  object,  its  image  at  the  sea  surface,  and  a  camera,  respec¬ 
tively.  The  points  R  and  S  are  respective  projection  of  P  and  Q  on  the  sea 
surface.  The  angles  of  incidence  and  reflection  are  designated  by  Q  *  and 
0  ,  respectively. 

Snell* s  law  yields  the  following  equations  for  X.  and  Y : 


X  -  Xo  —  W0  —  ----- — -  -  -tav»  0 

R$T 

V-  tyo  —  — Q~r rJr f-ttn  Q 

R  S 

in  which  the  distance  RS  is  given  bys 

=  -  Z  lavi  0  '  +  Z0  -ta 


and  there  is  the  following  relationship  between 

s  i  ^  0 


•;  i  v>  9' 


-  TL 


v\  Q 

Q  andQ 


(1) 

(2) 

(3) 

(it) 


in  which  n  indicates  the  refractive  coefficient  of  the  sea  water.  When  the 
position  of  the  camera  Q  is  fixed,  the  image  P*  can  be  determined  from 
equations  (1)  and  (2)  by  moving  the  point  P  on  the  object.  In  equations  (1) 
and  (2),  0  is  determined  as  a  function  of  P(x,y,z)  from  conditions  (3)  and  (1±). 

In  order  to  estimate  the  curvature  of  the  image,  let  us  consider  a  simple 
example  in  which  a  submerged  straight  bar  extends  with  an  angle  from  the  sea 
surface , (Figure  1|).  The  origin  of  coordinates  is  taken  at  the  intersection 
of  the  bar  with  the  sea  surface.  It  is  assumed  that  the  camera  is  located 
in  the  plane  perpendicular  to  the  bar,  and  its  coordinates  become  (0,L,0/L). 
When  the  angle  of  the  bar  with  the  sea  surface  is  taken the  coordinates 
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of  the  bar  are  given  by  (x,0,'$x).  Eliminating  x  from  (l)  and  (2)  under  the 
constraint  conditions  (3)  and.  (It),  we  have: 


xa+  C  L-Y 


(5) 


as  an  equation  for  the  curve  of  the  image.  At  first  sight  this  equation  does 
not  represent  a  straight  line  as  is  the  underwater  object.  However,  actual 
calculation  of  equation  (5>)  indicates  that  the  curve  expressed  by  (5>)  is  al¬ 
most  a  straight  line  for  nominal  conditions  of  aerial  photography  of  dye 
patches,  as  the  following  numerical  example  indicates.  Figure  £  shows  the 
computed  image  for  Q(  =  1,  0*1,  and  n  =  ,/2  ,  that  is,  when  the  angle  of 

the  camera  with  the  object  equals  lH?0  and  the  submerged  bar  makes  an  angle 
of  £*7°  with  the  sea  surface. 

When  it  is  assumed  that  the  curved  trail  of  dye  patches  is  caused  by  a 
vertically  differential  advection,  we  can  compute  distributions  of  vertical 
eddy  viscosity  of  the  mean  current,  if  the  curvature  of  the  trail  and  the 
velocity  and  shear  of  the  mean  flow  are  known  at  each  depth  (Ichiye,  1963). 
The  derivation  of  such  relation  is  briefly  discussed.  In  Figure  6  the  dis¬ 
tortion  of  a  dye  filament  due  to  the  differential  advection  is  schematically 
shown.  The  origin  of  coordinates  is  taken  at  the  sea  surface.  The  hatched 
areas  represent  the  dye  filament  at  the  initial  time,  and  after  time  t.  The 
coordinates  C  and  are  projection  of  the  dye  filament  on  the  sea  surface. 
The  components  of  the  mean  flow  are  expressed  by  u(x,y,z,t)  and  v(x,y,z,t). 

In  order  to  simplify  the  problem,  it  is  assumed  that  the  diffusion  of  the 
dye  filament  is  neglected.  Then  the  movement  of  the  filament  is  only  due 
to  the  advection  by  the  mean  flow,  and  there  are  the  following  relations: 

44  =  u  (>, '-y  z,t)  ;  y,  z,  t) 


(6) 
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The  radius  of  the  curvature  of  the  projected  curve  is  designated  by  R(s), 
in  which  ds  is  a  line  element  defined  by: 


(7) 


A  theorem  of  differential  geometry  yields: 


R  1  (A)  =  I  dzrt/d^ 


(8) 


-> , 


in  which  is  the  radial  vector  (£  ,  '?).  Since  it  is  more  convenient  to 
express  R(s)  as  a  function  of  the  argument  z,  we  have  the  following  equation: 


Riz)'1  - 


^  -  C 

/ 7 


(9) 


In  the  stationary  Ekman  flow,  the  components  of  the  mean  current 
satisfy  the  following  equations: 


(  k.  -■  d  5  i/  (  t 7  \  '  L  U.  i 


^  u "  A  j  iy 


(10) 


in  which  K(z)  is  the  vertical  eddy  viscosity  variable  with  depth  and  f  is 
Coriolis*  parameter.  The  values  of  ^  and  h  ,  obtained  by  integrating  (6), 
simply  equal  the  velocity  multiplied  with  time,  because  u  and  v  are  dependent 
only  on  z.  Substitution  of  these  values  of  and  into  (9)  yields: 


KW  =  f  R(z)  l(u*H  v‘)  1  (u‘ t  u  ')  V 


(11) 


This  equation  indicates  that  we  can  determine  the  vertical  eddy  viscosity  at 
any  depth  if  we  know  the  curvature  of  the  projection  of  dye  filament  R(z), 
the  mean  flow  u  and  v,  and  the  shear  uz  and  vz. 

Application  of  equation  (11)  to  the  dye  patches  in  the  previous  experi¬ 
ments  is  limited  because  of  uncertainty  of  the  depth  of  the  trailing  part. 
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However,  if  we  assume  that  the  eddy  viscosity  is  almost  uniform  in  the 
upper  ten  meters,  the  aerial  photographs  taken  at  1326  on  October  19,  1961 
(ENY,  GK)  can  be  used  for  estimation  of  the  eddy  viscosity.  In  this  example, 
the  velocities  u  and  v  have  an  order  of  magnitude  of  10  cm/sec,  the  shear  uz 
and  vz  of  ICT^sec-l,  t  is  10^  sec,  and  the  average  radius  of  curvature  esti¬ 
mated  from  the  photograph  is  about  10^  cm.  When  these  values  are  used, 
equation  (11)  yields  the  vertically  averaged  value  of  K(z)  equalling  about 
10^  cm^/sec,  which  is  quite  reasonable  compared  with  other  estimation  of  the 
vertical  eddy  viscosity  (Ichiye,  1962b). 

V.  STRIATION  OF  DYE  PATCHES 

Another  remarkable  feature  observed  in  many  experiments  is  striation 
(sometimes  split  into  several  blocks)  of  dye  patches  ('EM;  Ichiye,  1962a). 

Two  possible  causes  of  this  phenomenon  were  discussed  in  the  previous  paper 
(Ichiye,  1963).  One  is  a  hypothesis  on  mechanism  of  natural  streaks  proposed 
by  Welander  (1963).  The  other  is  the  roll- type  instability  in  rotating  tank 
experiments  by  Faller  (1963).  The  striation  of  dye  patches  has  more  regular 
and  wider  spacings  (an  order  of  100  m)  and  deeper  penetration  (a  few  tens  of 
meters)  than  natural  streaks  or  Langmuir  streaks.  Therefore,  WelanderTs 
model  is  hardly  applicable  to  the  present  phenomenon.  In  fact,  Reynolds* 
number  of  a  rotating  boundary  layer  ry  equals  about  300  for  a  number 

of  experiments,  in  which  v  10  cm/sec,  s  0.1  cm^/sec,  and^}__  =  I0“^sec“^. 
The  fact  that  this  value  is  far  above  the  critical  Reynolds*  number  15>0,  above 
which  the  roll-type  instability  starts  in  Faller* s  experiments,  seems  to  be 
favorable  to  the  cause  due  to  instability. 

The  third  possible  cause  of  striation  is  the  effect  of  waves  producing 
rolling-up  of  a  filament  of  dye,  as  discovered  by  Hama  (1962 )  by  numerical 
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integration  of  streaklines.  He  showed  that  the  streaklines  (filaments  of 
dye  injected  constantly  from  a  fixed  source)  become  rolled  up  into  a  number 
of  apparent  discrete  vortices  near  the  critical  layer  where  the  mean  shear 
flow,  given  by: 


a) 


has  zero  curvature  (y  e  o),  due  to  unamplifying  travelling  sinusoidal  waves 
superimposed  on  this  shear  flow.  It  is  noted  that  many  cases  of  striation 
were  observed  when  there  were  swells  and  winds  were  rather  high.  This  sug¬ 
gests  that  the  swells  and  the  wind  drift  may  provide  the  same  effect  as  the 
shear  flow  and  the  perturbation  superposed  on  it,  as  postulated  in  Hama!s 
model.  Also,  since  the  main  part  of  the  dye  patches  was  large  compared  with 
individual  striated  patterns,  the  dye  is  considered  to  be  fed  continuously 
to  the  streaklines.  It  would  be  an  interesting  problem  to  test  the  applica¬ 
bility  of  Hama*s  theory  to  striation  by  future  field  experiments. 
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